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ABSTRACT 

Non-blazar AGN have been recently established as a class of gamma-ray sources. M87, a nearby representative of 
this class, show fast TeV variability on timescales of a few days. We suggest a scenario of flare gamma-ray emission 
in non-blazar AGN based on a red giant interacting with the jet at the base. We solve the hydrodynamical equations 
that describe the evolution of the envelope of a red giant blown by the impact of the jet. If the red giant is at least 
slightly tidally disrupted by the supermassive black hole, enough stellar material will be blown by the jet, expanding 
quickly until a significant part of the jet is shocked. This process can render suitable conditions for energy dissipation 
and proton acceleration, which could explain the detected day-scale TeV flares from M87 via proton-proton collisions. 
Since the produced radiation would be unbeamed, such an events should be mostly detected from non-blazar AGN. 
They may be frequent phenomena, detectable in the GeV-TeV range even up to distances of ~ 1 Gpc for the most 
powerful jets. The counterparts at lower energies are expected to be not too bright. M87, and nearby non-blazar AGN 
in general, can be fast variable sources of gamma-rays through red giant /jet interactions. 

Subject headings: active galactic nuclei: jets - TeV photons: variability - stars: red giant 



1. INTRODUCTION 

Active galactic nuclei (AGN) are believed to be pow- 
ered by an accreting supermassive black hole (SMBH) 
in the center of a galaxy, a significant fraction of AGN 
show powerful jets, supersonic relativistic flows, on small 
(sub parsec) and larg e (multi hundred kpc) scales, (e.g. 
Begelman et al. 1984). These AGN are characterized by 



sion from one or many clouds from the BLR interacting 
with a hydrodynam ical jet recently has bin analyzed by 



Araudo et al. (20101. The radiation from the interaction 
between a massive star and a n AGN jet was studied in 



nonthermal emission extending from radio to high en 
ergy gamma-rays. This radiation comes from an accre- 
tion disc and from two relativistic jets that are launched 
close to the SMBH in two opposite directions. The emis- 
sion associated to the accretion process can be gener- 
ated by thermal plasma in the fo rm of an optically-thick 
disc under efficien t cooling (e.g., Shakura||l972 Shakura 
fc Sunyaev 1973), or as an optically- tmn corona (e.g., 
Bisnovatyi -Kogan fe Blinni kov| [T977[ [Liang fc Thomp~ 
son 1979). The emission from the jets is non thermal 



and ~comcs from a population of relativistic particles ac 
celerated for instance in strong shocks, although other 
scenarios are possible as well 
19981 



Bednarek fc Protheroe (1997). Namely they suggested 
that the jet interacts with stellar winds of massive stars, 
in their model they assume that the source of gamma- 
rays is moving with a relativistic speed, therefore the 
radiation is Doppler boosted. The main radiation mech- 
anism in this scenario is related to the development of the 
pair cascade in the field of the radiation of the massive 
star. 

In this work, we study the interaction of a red gi- 
ant (RG) star with the base of the jet in AGN and 
their observable consequences in gamma rays. We fo- 
cus here on the case of M87, a nearby non-blazar AGN 
that presents very high-energy r ecurrent activity with 
variability timescales of few days (lAharonian et al.|2006 



see, e.g., Schopper et al. 



Neronov fc Aharonian 2007 Rieger et al. 



Albert et a l. 2008; jAcciari et al:||2009[ |2010[ ). In the 



Rieger fc Aharonian||2008[ ) . This non-thermal ernission 
is thought to be produced through s ynchrotron and in - 



verse Compton (IC) processes (e.g. |Ghisellini et al.|1985 ) 
although hadr onic models have been also considered in 
the past (e.g.,|Manj iheimlll993l I Aharonian]|2000| |Mucke 
fc Protheroe]|20UTf pVharonian 2002^. 

'I'he existence of a stellar clustering in the central re- 
gions of AGN, possibly dow n to very s mall distances from 
the central SMBH (e.g. Penston 1988), implies that the 
interaction between a star and the jet should eventually 
happen. The gamma-ray production due to the inter- 
action between an obstacle and an AGN j et has been 
studie d in a number of works. For instance. Par fc Laor 
(1997) suggested the high-energy radiation produced by 
a beam of relativistic protons impacting with a cloud 
of the broad-line region (BLR). The gamma-ray eniis- 



tramework presented here, the jet impacts the RG en 
velope, already partially tidaly disrupted by the gravi- 
tational field of the central SMBH. The RG envelope is 
blown up, forming a cloud of gas accelerated and heated 
by jet pre ssure. The jet base is likely strongly magne 



tized (e.g., Komissarov et al. 2007 Barkov & Komissarov 



20081). The jet flow attected by the impact with the RG 



envelope can be a suitable region for particle accelera- 
tion, and a significant fraction of the involved magnetic 
and kinetic energy of the jet can be transferred to protons 
and electrons. Although electrons may not able to reach 
TeV emitting energies because of the expected large mag- 
netic fields, protons would not suffer from this constraint. 
These protons could reach the star blown material, and 
optically-thick proton-proton (jpp) interactions could lead 
to significant gamma-ray produc tion in the early stages 
of the cloud expansion. Unlike in [Bednarek fc Protheroe 
( 19971, we deal with solar-mass-type stars instead of the 
more rare high-mass stars, study the RG atmosphere-jet 



interaction, and follow the hydrodynamical evolution of 
the cloud. Finally, we do not introduce any beaming fac- 
tor to the radiation, since in our scenario most of the 
emission is produced when the cloud has not been sig- 
nificantly accelerated, Doppler boosting being therefore 
negligible. 

2. THE MODEL 

Main sequence stars are too compact to be signifi- 
cantly affected by tidal forces from the SMBH, unlike 
RGs, whose external layers are far less gravitationally 
bounded to the stellar core. Therefore, in the vicinity of 
a SMBH, the external la yers of an RG will suffer signifi- 
cant tidal disruption (se elKhokh lov et al.|1993b|a||Diener 



|Ayal et aL] r2000; Iva nov et aT 



odato 



etan mm 

et al.| |2OO90^ which can unbound from the stellar core a 
cloud with significant mass > 10'"' g. Therefore, if an RG 
penetrates into the innermost region of the jet, the RG 
envelope can be already weakly gravitationally attached 
to the star due to tidal disruption. In this situation, the 
external layers of the star can be lost due to jet ablation, 
which is unlikely in the case of undisrupted RGs (except 
for very powerful jets). 

The tidal forces are important when the distance be- 
tween the SMBH and the star is similar to or smaller 
than the tidal distance (zt) for a given RG radius (Rrg) 
and mass (Mrg), where: 



zt = R 



RG 



Mbh 
A/rg 



1/3 



and AfsH is the mass of the SMBH. Therefore, for a given 
RG-jet interaction distance to the SMBH z, the RG can 
lose the atmosphere layers beyond R*t- For the case of 
M87, with A^ RH = (6.4 ± 0.5) x 10^ AT© ( |Gebhardt fc" 
Thomas||2009), one obtains: 



-^RG ~ ^ 



\MbhJ 



1/3 



76 M^^^Rq 



(2) 

where Mrg© = ATrg/I Af©- Since a solar-mass RG, the 
most common one, can have up to few hundreds of Rq , 
a significant fraction of the star envelope can be carried 
away by the jet flow up to z < 10^^ cm. Note that 
evidence for the presence of a radio jet has been found 
from M87 within a distance of ^ 10^^ cm from the SMBH 



( Junor et al. 



'i'he M87 ' 



1999) 



'eV lightcurve obtained by | Aharonian et al 



(2006) shows several peaks, and each of these peaks in 



our model correspond to different RG-jet events. Note 
however that some nearby peaks may correspond to a 
complex disruption process, motivated for instance by a 
very disrupted and massive envelope, or by jet inhomo- 
geneities. Also, it cannot be discarded that a cluster of 
several RGs could also enter the jet. 

The time needed by the RG to cross the jet can- 
not be shorter than the typical M87 event duration, 
ie ^ 2 X 10^ s. It cannot be longer either, since then the 
event duration would also be longer if there is available 
RG matter for removal, as expected at z < zt. There- 
fore, tjc = te, and the interaction height can be derived 
taking the velocity of the RG orbiting the SMBH as the 



Keplerian velocity: 



GAL 



BH 



1/3 



-^Ql6g,_2/3^^ 



(3) 



where 0_i = 6'/0.1 is the jet semi-opening angle in ra- 
dians. An important para meter is the powe r of the jet 



U 



1 



10 



"'^ erg s ^ ( Owen et al.||2000 l, which we 



fix to Lj w 2 X 10 erg s~ . From Lj and the jet width, 
Zjcd, we can derive the jet energy fiux at the interaction 
height: 

(4) 



Fi 



ii 



TTZ.^t 



lO^'^erg cm ^ s ^ 



There are two regimes for the RG tidal disruption: un- 
der strong tidal interaction (i?RG > ^*t), the RG enve- 
lope suff ers an elongation along t he direction of motion of 
the star ( Khokhlov et al.|1993a ); under weak tidal inter- 
act ionJ^Kg^l^^lTlIHilenvelope is still roughly spheri- 
cal (Khokhlov et al. 1993b I. In both situations, the outer 



layers ot the star will be swept away by the jet, forming 
a cloud that will quickly heat up and expand. We study 
the time evolution of the cloud adopting a very simpli- 
fied hydrodynamical model for the cloud expansion. The 
heating of the cloud is caused by the propagation of shock 
waves, which are formed by the pressure exerted by the 
jet from below. Therefore, the cloud pressure is taken 
similar to the jet pressure (regardless it is of kinetic or 
magnetic nature): 



Pi 






{j-l)c 



(5) 



S.SxlO^^Mp/joCm, 



where c is the speed of light, and 7 the adiabatic in- 
dex (7 = 4/3). The cloud expands at its sound speed 
(cs), since the lateral and top external pressures are much 
smaller than the jet bottom one. When the cloud has 
significantly expanded, its pressure becomes smaller than 
the jet pressure from below. At that point new shocks 
leading to further cloud heating. We illustrate in Fig. [T] 
for the simplest case of weak disruption, how the spher- 
ical cloud evolves under the effect of jet pressure as seen 
in the plane perpendicular to the jet axis. 



2.1. Weak tidal interaction (spherical case) 

tions that characi 
;an be written as 1 

47rrj? 47rFjrj? 



The system of equations that characterizes the weak 
tidal interaction case can be written as follows: 



En = er 



3(7 -l)c 



drc 
dt 



crzc 
dt^ 



l{l-l)Ec 
cA€ ' 



1/2 



(6) 

(7) 
(8) 



where re and M^ are the cloud radius and mass, respec- 
tively. 
The solutions to Eqs. (6]|8) are: 



Vr{t) = 



TcO 



(l-iAcc)2 

2rco 



(9) 
(10) 




Fig. 1. — Sketch of the evolution within the jet of the cloud 
formed by the disrupted envelope of the RG. The plane of the 
image would be the jet section. 




Fig. 2. — Sketch of the proton acceleration and gamma-ray pro- 
duction processes. The plane of the image would be normal to the 
jet section. 

where Vco, assumed to be similar to R^.^, is the initial 
cloud radius, and ico the cloud characteristic expansion 
time: 



3cMr 



\T^lFjrcQ 



1/2 



5 X 10^ (Me28/^jM4r-c0.13)'/' 



(11) 

where Mc28 = A^c/lO^^g. Neglecting the initial cloud 
velocity in the z-direction, we obtain: 

m r,o { t Y 3/2-t/t^ 

Vz(t)^— 7Z -T-T^^ (13) 




Fig. 3. — Evolution of re (solid line) and z — Zj^ (dot-dashed 
lime) with time in the weak tidal interaction case. The parameter 
values are characteristic of M87: Lj = 2 X 10'*'* erg s~^, Mbh = 
6.4 X 10^ Mq, 9-1 = 0.5, Mrg = 1M0, 2jc Si 2.5 x 10^^ cm, and 
Mc Si 1.3 X 10^8 gr. 

where Zjc is the RG-jet penetration height. The evolution 
of the cloud radius is presented in Fig. [3] The adopted 
parameter values are: Lj = 2x lO^** erg s , Mbh = 6.4 x 
10^ Mq, e_i = 0.5, Mrg = IMq, Zjc « 2.5 x 10^^^ cm, 
Mc w 1.3 X 10^^ gr. Note that for times t < tec, z-Zjc < 
re < 9 Zjc and Vz <^ Cg <^ c. 

Making the cloud and jet pressures comparable, the 
energy transfer can be overestimated beyond a certain 
radius (^ct) and time (it) during the cloud evolution, in 
which energy balance is to be fulfilled: 



dt 



< TTrlF; . (14) 



Using Eqs. ([9| and (10), (14 1 permits to derive tt'- 



tt 



in 



1 



rco 



7 - 1 ic( 
Substituting Eq. (15) into Eq. (lo]), we obtain ret: 

''cO . -, r . . -,^^14,^^1/3 rn-l/3 



8 r,o 
7-1 tcoC 



2/3 



1.5 X io"m,^//8f;,4 



(15) 



(16) 



Thus, if Tc < r^t, the solutions presented in Eqs. (|6])- 
([8| are valid. After it, dE^/dt ~ Trr^i^j yields a slower 
increase of r^ with i, although this should be still a fast 
exponential growth. We consider i < it since we focus 
here on the case when the cloud is optically thick to pp 
collisions (see below), and at t > it the cloud density is 
already too low. 

2.1.1. Strong tidal interaction (elongated case) 

In the case of strong tidal interaction the RG atmo- 
sphere is stretched in the direction of motion of the star, 
and the expansion will be now cylindric. In such a case, 
-^RG ~ '^co ('''c is the cloud cylindrical radius) can be 
significantly smaller than the length of the disrupted at- 
mosphere, /c (Ayal ct al._2000J . The system of equations 
describing this case can be written as: 



7ic 



(1 - i/ice)3 



Er 



(7-l)c 



(17) 



drc 
~dt 






Me 
^ cMr 



1/2 



Substituting Eq. (17) into Eq. (18), we obtain 



rc{t) = reoe*/*^ 



(18) 
(19) 

(20) 



As in the weak case, r^o and ice are the initial radius and 
the expansion time of the cloud, where: 



with 



/ cMc 



V7r7fj/c 



1/2 



l/2p-l/2,-l/2 



Il4 



= lK;28^j,14 



c,14 



'day, (21) 



and F\ 



j,14 



(Fj/lO^'^ergcm-^s-i) 

If we neglect the initial velocity of the cloud in the z- 
direction, the distance covered by the cloud is as follows: 



Z{t) - Zjc ^ 

with a velocity 



2Fjlcrcotc 



cMr. 



1 t^cc^ 



tc 



vM = 



2Fjlcrcotc 
cMr 



(e*/*" - l) . 



(22) 



(23) 



As in the weak case, after substantial expansion equal- 
izing cloud and jet pressures overestimates the energy 
transfer from the jet to the elongated cloud, and the re- 
lation dEc/dt ^ "''^c-^j should be used. This phase is 
characterized by a slower, but still quite fast, power-law- 
like expansion rate. 

3. RADIATION 

Particles could be accelerated in the shocked jet region 
below the cloud. As noted in Sect.[TJ the jet is probably 
magnetically dominated at z < zt. Therefore, one can 
estimate the magnetic field in the jet as follows: 



Si 



8i 

cz 



2fl2 



2QQ L]'lz^i ez\G, 



(24) 



where ij,44 = ij/10''^erg s~^. The expected magnetic 
field in the shocked jet region should be also strong, prob- 
able of a similar strength to -Bj . Under such a magnetic 
field, one can estimate the acceleration timescale: 



E 



E, 



qBjC 



0.l^E2Br^s 



(25) 



where ^ is the acceleration efficiency parameter, q is the 
particle charge, E2 = E/W TeV, and Bj,2 = BjlQi^ G, 
the maximum energy of protons and electrons are 



E. 



p max 



gSj re « 10" Sj,2 re.i4r'/' TeV (26) 



and 



E,, 



qc 



^Ss-Bj 



lOB.^^^C^/^TeV, 



(27) 



respectively, where ag = 1.6 x 10^'^. Equation (26 1 
is obtained from limiting the proton acceleration by 
Bohm diffusion escape from the interaction region, of size 



Tc,i4 = (''c/10^'* cm), and Eq. (27l is obtained from lim- 
iting the electron acceleration through synchrotron cool- 
ing. Even taking a high ^ ~ 10 (for mildly relativistic 
shocks, as those of supernova explosions, ^ ~ lO"*), elec- 
tron energies will be too low to explain the HESS spec- 
trum o£M87 up to energies of few 10 TeV (,Aharonian| 



et al.| 2006 ) , whereas protons may be accelerated up to 
ultra-nigh energies. In addition, the expected -Bjet-values 
could easily suppress any IC component. We note that 
even for diffusion faster than Bohm, or under bigger ^- 
values, protons could still reach enough energy to explain 
observations. On the other hand, the cloud density can 
be high, making of pp interactions the best candidate for 
gamma-ray production in the RG-jet scenario, the char- 
acteristic cooling time for pp collisions being: 



IQi 



pp 



= l0^n's, 



(28) 



where rie. 10 = ric/lO cm^^ is the cloud density. We 
note that the high cloud density should not affect sig- 
nificantly the proton acceleration, which would occur in 
the far less dense jet shocked region. Nevertheless, pro- 
tons should penetrate in the acceleration process and, in 



the Blanford-Znajek scenario of jet f ormation (Blandford 



& Znajek 
formed on 



1977 Beskin et al. 



y by pairs at Zjc- 'I' 



1992) the jet is probably 



leretore, some cloud mate- 



rial should penetrate into the shocked jet medium, which 
can occur through Rayleigh- Tay lor instabilities (Chan^ 
|drasekhar||1961[ |Imshennik||1972[) . We present in TigT^ 
a sketch of the mixing, proton acceleration and gamma- 
ray production processes. We do not specify here the 
physics of particle acceleration, although this could take 
place by one or a combination of different mechanisms: 
magnetic reconnection right after the shock in the jet, 
shear acceleration due to the strong velocity gradients 
close to the contact discontinuity, or some sort of stochas- 
tic acceleration due to magnetic turbulence downstream 
the jet shock. Regarding other proton radiation mech- 
anisms, proton synchrotron will not be efficient in our 



case, with tp sync « 5 x 10 " B-2 s^ tpp. Also photome- 

5 X 

1^ 



son production can be also neglected, since tp^ 
10^ L4ir~l^£~^ya > tpp, where L41 = {L/lO^^eTg s 
and ekov = (e/1 keV) would be the luminosity produced 
in the region and the ambient photon energy (e.g. ther- 
mal X-rays; see below), respectively, photomeson pro- 
duction with keV ambient photons would require protons 
with energies above ^ 100 TeV 

Hereafter, we will treat the generation of protons with 
energies > 100 GeV phenomenologically, assuming that a 
fraction 77 of the total eclipsed jet luminosity is converted 
to relativistic protons: Lp = rjTrr^Fj. We also assume 
that these protons can effectively reach the cloud, where 
they suffer pp interactions that lead to 7r°-meson produc- 
tion, although some of them could escape surrounding 
the cloud. Once in the cloud, protons can be effectively 
trapped during the RG-jet interaction time for magnetic 
fields as low as few G, since: 



Br = 



tEc 
iqr^ 



O.OSis-Bare^i^G, 



(29) 



where t^ = t/lO^ s is the time inside the cloud. Note 
that Eq. ( 29 ) has been derived assuming Bohm diffusion. 



but faster diffusion regimes would still allow to keep pro- 
tons trapped in the cloud. 

The typical fraction of the proton energy transferred 
per collision to the leading gamma rays is E^ = Q.YIEp 



( Kelner et al. 2006 ) in the optically thin case, and around 
twice that value for optically thick media neglecting 
gamma-rays from other secondary particles. Therefore, 
we can characterize the proton-gamma ray energy trans- 
fer by 



X = E^/Ep = 0.17 [2 - exp{-t/tpp)] . 



(30) 



Two phases of the cloud expansion can be distin- 
guished: the radiatively efficient regime, i.e. with x ~ 
0.34 or i > ipp, and the radiatively inefffcient regime, 
with X = 0-17 OT t < tpp. Thus, from the simplifications 
above, the gamma-ray luminosity in the pp optically- 
thick case can be written as: 



0.3477 7rr^Fj, 



(31) 



where is seen that L^ oc r^. In the pp optically thin case, 
only a fraction t/tpp of Lp is lost through pp collisions, 
and L^ ex r'^^. The general expression for the gamma-ray 
luminosity during a RG-jet interaction event becomes: 



^7 ~ ^VX r[ 



Fj (1 - eM~t/tpp)) 



(32) 



Given the fast expansion of the cloud, either in the spher- 
ical or the elongated case, one can expect a sharp spike 
in the light curve. 

Secondary electrons and positrons (e^), injected hy pp 
collisions with an energy rate ~ L^, could emit most of 
their energy through synchrotron radiation. Given the 
moderate energy budget, the radio, optical and X-ray 
fluxes would be below the observed values in the region of 
interest. However, at later times, conditions may change 
becoming more suitable for radio emission. Thus, it can- 
not be discarded that RG-jet interactions could eventu- 
ally have a low-energy faint counterpart produced by sec- 
ondary or thermal e , or even by a primary population 
of accelerated electrons. 

The gamma-ray pp light curve for M87 is presented 
in Fig. HI in which the maximum is reached at tpcak ~ 
4 X 10^ s, with a width of ^ 1 — 2 days. A value for ry 
of 0.1 has been adopted. We recall that our approach 
is valid for i < it ~ 4.7 x 10^ s (see Sect. [2|, hence 
the adopted cloud evolution model describes the gamma- 
ray peak properly. We show the gamma-ray lightcurve 
for the weak tidal disruption case as the most conser- 
vative scenario. In the strong tidal disruption case, the 
lightcurve would be similar but the gamma-ray maxi- 
mum would be even higher. 

3.1. Optimal radiation case 

Around the gamma-ray maximum, at i ^ ipeak, and 
because of ipeak ^ tpp, one can measure the cloud density 
and radius Ucp and rep. From 



3 Me 



"-^ An mp r3p 



plus the variability time: 



6v ^^ tp 



10 



15 



(33) 



(34) 




Fig. 4. — Gamma-ray pp lightcurve for the weak tidal disruption 
case. The same parameter values as in Fig. Is] have been adopted. 

one can determine 



_ 3Mc tv 
"P ~ ^ 47r 1015 m. 



1/3 



(35) 



We can characterize the variability time iv, which 
would correspond to the characteristic duration of the 
gamma-ray peak, as follows: 



iv = 2- ^^ 



t-cc V ^ ^pcak / t^cc / 



dL^/dt 



(36) 



The condition r ~ rep, and Eqs. (|9|), (35 1 and (36), allow 
the derivation of the following expression: 



48A/c 
TT 10^5 m^ 



1/3 



7/3 



f" r 



(37) 



Then, from Eq. (37), and assuming ico — ^jc, (1 
ipcak/tce) < 1 and rco = R^q, one obtains: 



''CO. '^^ 



^^3/20^3/2001/10^21 



/20 



50 X 10l/47r3/20c3/20G'l/5JV^^A0^3/20^^/10 



M, 



m^Mnotl 



1/5 



750'^Trc'^G^M^^mpt 



Ml,^G^Ll'^^/Hf'^ 



1/15 



3.16 X 10227r3/2m^/^c3/2M, 



(38) 



(39) 



(40) 



RGC 



Finally, substituting ^, (|4|, (|35| and (|39| into pTj ), 
the expression for the gamma-ray luminosity around the 
lightcurve peak can be derived: 

/ ,,r X 4/15 / rSrSi \ 1/5 



'' 1097r3/5 \M^^J I c3G4^35)8 



«8 X lO''.l-,L]%tl[lM'^^'M-^l''ezre.g s^ 

where Mbh,9 = (Mbh/IO*^ M©), 77_i = 77/0. 1 and iv,5 = 
(iv/lO^s). Adopting typical parameter values for M87, 
Lj = 2x104-* erg s-i,Mbh = 6.4 x109Mo,6I_i = 1, iv = 
2 X 10^ days, Mrg = IM©, Zjc « 3.6 x 10*^ cm, M^ « 



6 



1.4x 10 ° gr, and ?y_i = 1, one gets Lj w 4x 10 ' e rg s 



in good agreement with observations ( Aharonian et al 



2006| [Albert et ai:i[2008| |Acciari et aL][lT09| |20l0| ). We 



remark that the tv-value would be in agreement with the 
observed event durations. 

3.2. Thermal radiation of the cloud and self-j"f 
absorption 

In the case of M87, near the peak of the very high- 
energy (VHE) radiation, i.e. Ucp ~ 10^° cm^'^ and 
10^^ cm, the cloud is optically thin to the radi- 



' cp 



ation produced by its own shocked plasma: 

Te-y = ^cp ricp (Tt ~ 0.6 < 1, 



(42) 



where ctt = 6.65 x 10^^^ cm^^ is the Thompson cross 
section. At the temperature of the shocked cloud, Tc ^ 
10^*^ K, the timescales for Coulombian thermalization 
through p — p and e — e scattering are te_£ 



^3/2 1 



tp-f 



1000 s (tep « 10^ T^^ Hcio ~ 10 s for p - e 
scattering). Therefore, the shocked cloud is thermalized. 
The main channel of thermal radiation is free-free emis- 
sion, with a photon mean energy < e > ^-^ kT c ^ 1 MeV 
and total lumin osity ( Berestetskii et al^|1971 Kaplan & 
PikerNer|[T979l ): 



Lx = 2.l X lO^^^T^/^n^V; 



cp 



10^1 erg s" 



(43) 



where Fcp = 47rr^ /3. The concentration of thermal pho- 



tons can be estimated as 



nx- 



Att ckTcr^, 



10^ 



(44) 



cp 



yielding an optical depth for photon-photon absorption 



at the energy of the strongest attenu ation (^ 



'1< 



0.2nxrcpCTT ~ 10'^ < 1 (|Aharonian||2004[ ) , 
being much smaller at 1 TeV. The tree-tree radiation 
should not show any thermal lines, presenting a very hard 
non-thermal X-ray spectrum. 

Fo r v ery powerful jets the condition presented in 
Eq. (42 1 is not fulfilled, the shocked plasma is radiation 
dominated and cooler, and Eqs. (43) and (44 1 do not 



apply. The cloud is then optically thick, with the radi- 
ation being a black body with mean energy of photons 



< e >w 3fcTb 



104/44 ^^BH,9 



2/15 ,-21/60 



t. 



77 



eV. The optical 

26/15 



v,5 



depth for gamma-rays is r- 
with the radiation being suppressed for energies -Eth ^ 
TOg c^/3 A; Tb ^ 50 GeV, where nie is the electron mass. 
Photon-photon absorption creates pairs with energies 
> Sth that cool down through synchrotron emission with 
spectral energy distribution eF^ ex e""^/^ below 10 keV, 
with the higher energy part of the spectrum softer, reach- 
ing MeV-GeV energies. Under t^j > 1 and reasonable 
magnetic fields, the synchrotron luminosity will be simi- 
lar to the absorbed gamma-ray luminosity. 

The X-ray flare detected from M8 7 almost simultane- 
ously with the VHE flare (see, e.g., Acciari et al. 20091 



may have been also produced at the KG-jet interaction 
This X-ray flare could be of synchrotron nature, with 
possible contributions from a primary electron compo- 
nent, secondary e^ from pp and photon-photon inter- 
actions, and thermal free-free radiation. Regardless the 



origin, the observed X-ray emission could have a coun- 
terpart at lower energies. If it came from the region 
of gamma-ray production, the spectrum should be quite 
hard to avoid too many optical photons that would lead 
otherwise to significant TeV photon absorption. Optical 
observations simultaneous with a gamma-ray flare could 
clarify this point. 

4. DISCUSSION AND CONCLUSIONS 

The total jet l um inosity can be inferred from observa- 
tions using Eq. (41): 



Lj = 8 X 10 



44 r5/8 )\#2/3 



^"41^, 



BH.9'' 



-5/8,- 



l/8,r-l/6 



'M, 



RG 



erg s 



(45) 
This formula weakly depends on the observables, being 
almost insensitive to Mrq, on the other hand hard to es- 
timate. This provides therefore a quite robust estimate 
of the jet luminosity with 77 as most unknown parame- 
ter. Actually, if ij were known, then rj could be also 
estimated. 

For the most powerful jets, L-^ would be limited by 
the jet size becoming L^ — xv^i- Taking for in- 



stance Lj '• 
w 2 X 10^5 ' 
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erg s 



Ly could be as high as 



77_i erg s~ . An improvement of a factor 
of several in the VHE sensitivity (e.g. through the forth- 
coming the Cherenkov Telescope Array -CTA-) would 
test our gamma-ray predictions for the whole RG-jet in- 
teraction process, including the early cloud expansion 
phase, allowing for a detailed study of the involved (mag- 
neto)hydrodynamics, particle acceleration, and radia- 
tion. 

We remark that, if a detectable gamma- ray flare with 
a duration of few days were to be produced in M87, in 
particular through pp interactions, the cloud should have 
a mass of ~ 10^* g. Such a massive cloud cannot ac- 
quire a large speed in the jet direction at the times when 
pp collisions are an efficient gamma-ray emitting mech- 
anism, and therefore the emission will not suffer signif- 
icant Doppler boosting. In the case of a lighter cloud, 
large Lorentz factors can be achieved, but then pp in- 
teractions will be inefficient producing gamma-rays, the 
probability to detect a flare lower due to beaming, and 
the duration of the event shorter than observed because 
of faster expansion and beaming. 

Coming back to the question of cloud mass, we note 
that to extract a cloud with a mass > lO^'^ g, a more 
powerful jet than in M87, for similar j et-R G interaction 
conditions, would be required (see Eq. p9[ )). 

An important question is whether there are enough 
RGs in M87 at the relevant jet scales. The model pre- 
sented here would require few interactions per year to 
explain the observations in M87. Since the typical dura- 
tion of the RG-jet interaction is of about 3-4 days, the 
RG fiUing factor should be T - 4/365 « lO'^. With a 
jet volume at the relevant scales of ~ ttO'^z^/S, the den- 
sity of RGs in the region should be ^ T/V ^ 2x 10^ pc^'^ 
for M87. Unfortunately, no direct information is avail- 
able on the density of stars in the vicinity of the SMBH in 
M87. The stellar mass in a sphere with a radius of 80 pc 
is est imated in 2 x IO^Mq (e.g. [Gebhardt fc Thomas 



2009| , and these observational data should be extrapo- 
lated four orders of magnitude down to ~ 0.01 pc. Thus, 
depending on the assumed extrapolation law, the num- 



ber of RGs in the vicinity of the SMBH may or may not 
be enough. It is worth noting that a dense stellar cluster 
near the SMBH could be behind the broad-line region 
in AGN as produced by the blown-up atmosphere of red 
dwarfs, which would imply the presence of numerous RGs 
in the center of AGN (Penston^l988} . In addition, stud- 
ies of the possible ste llar density profiles in the vicinity of 
the SMBH in A GN ( |Bisnovatyi-Kogan et aLl|1982[ |Mur ' 



phy et al.|1991 1 show that densities liJie the required one 
(^ 2 X 10 pc^) could be achieved. The observation of 
VHE flares could be already an indication that enough 
RGs are present near the SMBH in M87. 

Interestingly, RG/jet interactions are expected to be 
transient phenomena. At higher jet heights, although 
many RGs could be simultaneously present in the jet 
rendering rather continuous emission, the much more di- 
luted jet would not remove a significant amount of ma- 
terial from the star and the effective cross section of the 
interaction would be just i?RG, yielding a low energy bud- 
get for such a multiple interaction events. 

The scenario presented here, adopted to explain the 
day VHE flares observed from M87, could also be rele- 
vant in other non-blazar AGN. For blazar sources the 
beamed emission would overcome the RG-jet interac- 
tion, expected to be weakly beamed due to moderate v^- 
values. For instance, the closest AGN, the radio galaxy 
Gen A, at -- 3.8 Mpc distance ( |Rejkuba| [20041 ), could 
also show detectable flare like emission. At p resent, per 



sistent fain t VHE emission has been detected (jAharonian 
et al.|2009[ ) with L^ = 2.6 x lO^^ 



erg s 



Accounting for 
the black hole mass of this AGN, Mbh = 5.5 x lO'^M©, 
taking the observed VHE luminosity as a reference, and 
assuming t 
a jet luminosity ij 
est value. Therefore, it cannot be excluded that RG-jet 



1 day, one derives implementing Eq. 45 
1.2 X lO"'^ erg s~^, a rather mo3~ 



interactions may contribute to the VHE radiation de- 
tected from Gen A, or that transient activity due to RG- 
jet interactions may be observed from this source. An- 
other cas e, the radio galaxy NGG1275, at a distance of 
73 Mpc fHick en et al . |2009), shows variable behaviour 



in GeV ( ,Abdo et al. 
about 2 X lO*"* erg s 



09). The GeV luminosity is of 



Using Eq. ( 45 1 , we can esti- 
10^^ 



mate the power of the jet as 5 x lO'^erg s adopt- 
ing an Mbh = 10*^ Mq. In the case of NGC 1275 the 
shocked cloud would be optically thick at the luminos- 
ity peak, implying significant attenuation of the TeV 
emission through photon-photon absorption with a cutoff 
around 50 GeV. 

At farther distances, the strong jet luminosity dcpcn 
deuce L^ oc 
L 
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L^-^ implies that FR II sources with say 
^ may be still detectable up to distances 



-J - -" erg s 
of ^ 0.5 Gpc (i nternal absorption shou ld be included in 



these cases; see Aharonian et al. 2008). Also, the lumi 



nosity in the range 0.1-100 GeV would also be significant 
unless there is a strong low-energy cutoff in the proton 
spectrum. Therefore, Fermi may detect day-long GcV 
flares originated due to RG-jet interactions from FR II 
galaxies up to distances of few 100 Mpc. Summarizing, 
GeV and TeV instrumentation can potentially detect a 
number of RG-jet interactions per year taking place in 
nearby FR II and very nearby FR I galaxies, with the 
most powerful events being detectable up to 1 Gpc. 
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